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The boundary-value problem for Laplace-type operators acting on smooth sec- 
tions of a vector bundle over a compact Riemannian manifold with generalized local 
p boundary conditions including both normal and tangential derivatives is studied. 

1^ ■ The condition of strong ellipticity of this boundary-value problem is formulated. 

The resolvent kernel and the heat kernel in the leading approximation are explic- 
itly constructed. As a result, the previous work in the literature on heat-kernel 
^ ' asymptotics is shown to be a particular case of a more general structure. For 

d • a bosonic gauge theory on a compact Riemannian manifold with smooth bound- 

ary, the problem is studied of obtaining a gauge-field operator of Laplace type, 
jointly with local and gauge-invariant boundary conditions, which should lead to 
a strongly elliptic boundary-value problem. The scheme is extended to fermionic 
gauge theories by means of local and gauge-invariant projectors. After deriving 
a general condition for the validity of strong ellipticity for gauge theories, it is 
proved that for Euclidean Yang-Mills theory and Rarita-Schwinger fields all the 
above conditions can be satisfied. For Euclidean quantum gravity, however, this 
property no longer holds, i.e. the corresponding boundary-value problem is not 
strongly elliptic. Some non-standard local formulae for the leading asymptotics of 
the heat-kernel diagonal are also obtained. It is shown that, due to the absence of 
strong ellipticity, the heat-kernel diagonal is non-integrable near the boundary. 
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1 Introduction 

The consideration of boundary conditions in the formulation of quantum field theories is 
crucial for at least two reasons: 

(i) Boundary effects are necessary to obtain a complete prescription for the quantization, 
unless one studies the idealized case where no bounding surfaces occur. In particular, the 
functional-integral approach relies heavily on a careful assignment of boundary data. 

(ii) One may wonder whether the symmetries of the theory in the absence of boundaries 
are preserved by their inclusion. Moreover, one would like to know whether the applica- 
tion of a commonly accepted physical principle (e.g. the invariance under infinitesimal 
diffeomorphisms on metric perturbations, in the case of gravitation) is sufficient to de- 
termine completely the desired form of the boundary conditions, upon combination with 
other well known mathematical properties. 

The investigations carried out in our paper represent the attempt to solve these prob- 
lems by using the advanced tools of analysis and geometry, with application to the oper- 
ators of Laplace or Dirac type. Indeed, the differential operators of Laplace type are well 
known to play a crucial role in mathematical physics. By choosing a suitable gauge it is 
almost always possible to reduce the problem of evaluating the Green functions and the 
effective action in quantum field theory to a calculation of Green functions (and hence 
the resolvent) and functional determinants (or the ^-functions) of Laplace-type operators. 
These objects are well defined, strictly speaking, only for self-adjoint elliptic operators. 
Thus, on manifolds with boundary, one has to impose some boundary conditions which 
guarantee the self-adjointness and the ellipticity of the Laplace-type operator. The sim- 
plest choice are the Dirichlet and Neumann boundary conditions when the fields or their 
normal derivatives are set to zero on the boundary. A slight modification of the Neumann 
boundary conditions are the so-called Robin ones, when one sets to zero at the bound- 
ary a linear combination of the values of the fields and their normal derivatives. An even 
more general scheme corresponds to a mixed situation when some field components satisfy 
Dirichlet conditions, and the remaining ones are subject to Robin boundary conditions 

However, this is not the most general scheme, and one can define some generalized 
boundary conditions which are still local but include both normal and tangential deriva- 
tives of the fields [Q. For example, by using linear covariant gauges in quantum gravity, 
one is led to impose boundary conditions that involve the tangential derivatives as well, 
to ensure that the whole set of boundary conditions on metric perturbations is invari- 
ant under infinitesimal diffeomorphisms 0, |^. The same boundary conditions may be 
derived by constructing a BRST charge and requiring BRST invariance of the boundary 
conditions 0. 

In this paper we are going to study the generalized boundary-value problem for 
Laplace- and Dirac-type operators, the latter being relevant for the analysis of fermionic 
models. Such a boundary-value problem, however, is not automatically elliptic. There- 
fore, we find first an explicit criterion of ellipticity (section 2.3). Then we construct the 
resolvent kernel and the heat kernel in the leading approximation (section 3). The appli- 
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cation of this formalism (section 4) proves that the generahzed boundary- value problem is 
strongly elliptic for Euclidean Yang-Mills theory (section 5) and Rarita-Schwinger fields 
(section 6), but not for Euclidean quantum gravity (section 7). The possible implications 
are discussed in section 8. 

2 Laplace-type and Dirac-type operators 

Let (M, g) be a smooth compact Riemannian manifold of dimension m with smooth 
boundary, say dM. Let g be the positive-definite Riemannian metric on M and g be the 
induced metric on dM. Let ^ be a (smooth) vector bundle over the manifold M and 
C°°(y,M) be the space of smooth sections of the bundle V. Let V* be the dual vector 
bundle and E : V ^ V* he a Hermitian non-degenerate metric, E"^ = E, that determines 
the Hermitian fibre scalar product in V. Using the invariant Riemannian volume element 
dvol (x) on M one defines a natural inner product (, ) in C°°{V,M), and the Hilbert 
space L^(V,M) as the completion of C°^(V,M) in this norm. Further, let be the 
connection on the vector bundle V compatible with the metric E, tr g = g ^ 1 he the 
contraction of sections of the bundle T*M ^T*M ®V with the metric on the cotangent 
bundle T*M, and Q be a smooth endomorphism of the bundle V, i.e. Q G End(\^), 
satisfying the condition 

Q = E-^Q^E = Q. (2.1) 

Hereafter we call such endomorphisms self-adjoint. We also use the notation A = E^^A^E 
for any endomorphisms or operators. Then a Laplace-type operator, or generalized Lapla- 
cian, 

F : C^iy, M) C^iy, M), (2.2) 
is a second-order differential operator defined by 

F = -tr,V™^V^ + Q, (2.3) 

where 

yT*A/®y ^ yT*M ^ 1 + 1 ^ v^, (2.4) 

and V^**''^ is the Levi-Civita connection on M. 

Let 1/ be a Clifford bundle and F : T*M End(V), F(0 = F^^^^, be the Clifford 
map satisfying 

F(ei)F(6) + r(6)F(ei) = 2^7(^1, 6)/, (2.5) 

for all ^1,^2 € T*M. Hereafter / denotes the identity endomorphism of the bundle V. 
Let be the Clifford connection compatible with the Clifford map. Then a Dirac-type 
operator 

D : C°^(V,M) ^ C°°(V,M), (2.6) 
is a first-order differential operator defined by [|l], 

D = i{TV^ + S), (2.7) 

with S G End {V). Of course, the square of the Dirac operator is a Laplace-type operator. 
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2.1 Geometry of boundary operators 
2.1.1 Laplace-type operator 

For a Laplace-type operator we define the boundary data by 

M^-i^). (2.) 

where 

ipoiip) = (p\dM, i^iiv) = Vatv^I^m (2.9) 

are the restrictions of the sections ip G C°^{V,M) and their normal derivatives, to the 
boundary (hereafter, is the inward-pointing unit normal vector field to the boundary). 
To make the operator F symmetric and elliptic (see sections 2.2 and 2.3), one has to 
impose some conditions on the boundary data ipp{ip). 

In general, a (i-graded vector bundle is a vector bundle jointly with a fixed decompo- 
sition into d sub-bundles In our problem, let the vector bundle Wp over dM be the 
bundle of the boundary data. Wp consists of two copies of the restriction of V to dM 
and inherits a natural grading 

Wf = Wo®Wi, (2.10) 

where Wj represents normal derivatives of order j, and, therefore, 

dimW^F = 2dim\/. (2.11) 

The bundles Wo and Wi have the same structure, and hence in the following they will be 
often identified. Let W'p = Wq © W[ be an auxiliary graded vector bundle over dM such 
that 

dim PF^ = dim K (2.12) 

Since the dimension of the bundle Wp is twice as big as the dimension of Wp, dimWp = 
2 dimiy^, it is convenient to identify the bundle Wp with a sub-bundle of Wp by means 
of a projection 

Pp: Wp^W'p, PI = Pf- (2.13) 

In other words, sections of the bundle W'p have the form PpX-, with x being a section of 
the bundle Wp. Note that the rank of the projector Pp is equal to the dimension of the 
bundle V , 

rank Pp = rank (Iw - Pf) = dim V, (2. 14) 

where Iw is the identity endomorphism of the bundle Wp. 

Let Bp : C^iWp, dM) C°°{Wj,, dM) be a tangential differential operator on dM. 
The boundary conditions then read 



Bpi:p{^) = 0. 



(2.15) 
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The boundary operator Bp is not arbitrary but should satisfy some conditions to make 
the operator F self-adjoint and elliptic. These conditions are formulated in the subsects. 
2.2 and 2.3. The boundary operator Bp can be then presented in the form 

Bp = PpL, (2.16) 

where L is a non-singular tangential operator 

L : C°°{WF,dM) C^iWp^dM), (2.17) 

meaning that there is a well defined inverse operator L~^. In other words, the boundary 
conditions mean 

^f{^) = KfX with arbitrary x e C^iWp, dM), (2.18) 

where 

Kp = L-\lw - Pf) ■■ C^iWF^dM) ^ C^iWF.dM). (2.19) 
Now let n be a self-adjoint projector acting on Wq and Wi, 

U: Wo^Wo, U: Wi^Wi (2.20) 

= n, 11 = E-^u^E = n. (2.21) 

In our analysis, we will consider the projector Pp and the operator L of the form 

p^ = n©(j-n) = (U (j°n))' ^^-^^^ 



/ 

A / 



(2.23) 



where A is a tangential differential operator on dM, A : C°°(Wo,dM) — > C°°{Wo,dM), 
satisfying the conditions 

nA = An = 0. (2.24) 
Of course, rank Pp = rank 11 + rank (/ — 11) = dim V as needed. Hence we obtain 

B.= ('! V Kp = I^-Bp=('-^-'' (2.25) 



,A iw-'nj' ^ ' V -A n. 

Note that Bp and Kp are complementary projectors 

Bl = Bp, Kl = Kp, BpKp = KpBp = 0. (2.26) 
Moreover, by virtue of the property ( p.l4|) , both the operators Bp and Kp do not vanish. 
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2.1.2 Dirac-type operator 

For a Dirac-type operator (see (|2.7| )) the normal derivatives are not included in the bound- 
ary data ipoif), and such boundary data consist only of the restriction ipoi'^) = ^0(9^) 
(see ( p.9D ) of sections ip G C°°{V,M) to the boundary. Therefore, the bundle of the 
boundary data Wd is just the restriction Wq of the bundle V and, similarly, the auxiliary 
vector bundle W'^y = Wq. The dimensions of these bundles are 

dim Wd = 2 dim W'^ = dim V. (2.27) 



This should be compared with ( p.ll| ) and ( p.l2| ). As above, W'j^ can be identified with a 



sub-bundle of Wd by means of a projection 

Pd: Wd^W'o, (2.28) 

but now the dimension of the projector is equal to half the dimension of the bundle 
V: 

rank = -dim y. (2.29) 
The boundary conditions for the operator D read 

BdMv) = 0, (2.30) 

where now the operator Bd is not a tangential differential operator but just a projector 
Pd: 

Bd = Pd, Kd = I- Pd- (2.31) 

It is crucial to use projectors to obtain a well posed boundary-value problem of local 
nature for Dirac-type operators, since any attempt to fix the whole fermionic field at the 
boundary would lead to an over-determined problem. 

2.2 Symmetry 

2.2.1 Laplace-type operator 

The Laplace-type operator F is formally self-adjoint, which means that it is symmetric 
on any section of the bundle V with compact support disjoint from the boundary dM, 
i.e. 

If{^i, ^2) = {Fipi, ip2) - (^1, Fip2) = 0, for v?i, e C^{V, M). (2.32) 

However, a formally self-adjoint operator is not necessarily self-adjoint. It is essentially 
self-adjoint if its closure is self-adjoint. This implies that the operator F is such that: i) 
it is symmetric on any smooth section satisfying the boundary conditions ( ^.15 ), i.e. 



/f(¥^i, ^2) = 0, for y^i, ¥^2 e C^{V, M) : PfV(</'i) = ^^^(¥^2) = 0, (2.33) 

and ii) there exists a unique self-adjoint extension of F . The property ii) can be proved 
by studying deficiency indices 0], but is not the object of our investigation. 
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The antisymmetric bilinear form Ip{ipi,ip2) depends on the boundary data. Integrat- 
ing by parts, it is not difficult to obtain 

If{Vi, V2) =< '^{'fi), JpipM >, (2.34) 

where 

Jf={^j o)- (2.35) 

Here <, > denotes the inner product in C°°(iy, (9M) determined by the restriction of 
the fibre metric to the boundary. 

To ensure the symmetry, we have to fix the boundary operators i^i? so as to make this 
form identically zero. The condition for that reads 

KfJfKf = 0. (2.36) 



By using the general form ( p.25 ) of the operator Kp, we find herefrom that this is equiv- 



alent to the condition of symmetry of the operator A 

< Aipi, >=< <^i, A<^2 >, for ^1, v?2 G C°°{Wo, dM), (2.37) 
Thus, we have proven the following result: 



Theorem 1 The Laplace-type operator F ( \2.3[ ) endowed with the the boundary conditions 
( ^.I3i ), with the boundary operator B given by ( ^.23i ), and any symmetric operator A 



satisfying the conditions ( \2. 24\) and ( \2. 3'1\ ), is symmetric. 



2.2.2 Dirac-type operator 

In the case of Dirac-type operators there are two essentially different cases. The point is 
that, in general, there exist two different representations of the Clifford algebra satisfying 

r(0 = -sm, e = ±1. (2.38) 

In odd dimension m there is only one possibility, e = —1, corresponding to self-adjoint 
Dirac matrices, whereas for even dimension m they can be either self-adjoint or anti-self- 
adjoint. By requiring the endomorphism 5* to satisfy the condition 

S = eS, (2.39) 

we see that the Dirac-type operator (|2.7| ) is formally self-adjoint for e = — 1 and anti-self- 
adjoint for e = 1: 

Id{vu ip2) = {D^i, ^2) + e(<^i, D^2) = 0, ipi, ^2 e C^{V, M). (2.40) 

In complete analogy with the above, we easily find the condition for the Dirac-type oper- 
ator to be (anti)-symmetric, 

Id{^i, ^2) =< ipoM, Jo-ipoM >= 0, (2.41) 
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Jd = ieT{N) = ieVN^, (2.42) 
for any fi,f2 € C°°(V,M) satisfying the boundary conditions (|2.3CI|) . This leads to a 



condition on the boundary operator, 

KdJdKd = 0, (2.43) 
wherefrom, by using (p.31|) , we obtain a condition for the projector Pd, 

{I-Pd)T{N){I-Pd)=0. (2.44) 
Hence we get a sufficient condition on the boundary projector, 

PD = T{N){I-PD)T{Ny\ (2.45) 
This means that Pd can be expressed in the form 

P^ = i(J + r^), (2.46) 

where r] satisfies the conditions 

r]'^ = J, f]T{N) + T{N)7] = 0. (2.47) 

Thus, there are two cases: i) r] is anti-self-adjoint and commutes with r(A^), ii) i] is 
self-adjoint and anti-commutes with r(A^). This leads to two particular solutions, 

7] = ±r(A^), for e = 1, and even m, (2.48) 

and ^ 

rj = - — :T(u), for e = —1 and any m, (2.49) 

\u\ 

with some cotangent vector u G T*dM on the boundary. 

In even dimension m there exists another very simple solution, 

1] = ±C, for even m, e = ±1, (2.50) 

where C is the chirality operator defined with the help of an orthonormal basis e" in T*M, 

C = z'"/2r(ei)...r(e'"). (2.51) 

One easily finds 

C2 = /, C = C, (2.52) 
CT(0 + mC = 0, (2.53) 

for any ^ G T*M, so that the conditions ( p.47| ) are satisfied. 
Thus, we have 
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Theorem 2 The Dirac-type operator D = ^(rV + S) with T and S satisfying the con- 
ditions T = —eT, S = eS, endowed with the boundary conditions PdvIsm = 0, is anti- 
symmetric for e = —1 and symmetric for e = 1 provided that the boundary projector 
Pd satisfies the condition (/ — Pd)^{N){I — Pd) = 0. Admissible boundary projectors 
satisfying this criterion are: for e = 1 



Pn = \[I±TiN)], (2.54) 



and, for e = —1, 

PdM 



(2.55) 



1 

2 

where u G T*dM. In the case of an even- dimensional manifold M, the boundary projector 



\u\ 



Pn = liI±C), (2.56) 
C being the chirality operator, is also admissible. 

2.3 Strong ellipticity 
2.3.1 Laplace- type operator 

Now we are going to study the elhpticity of the boundary-value problem defined by the 
boundary operator ( p.25|) . First of all we fix the notation. By using the inward geodesic 
flow, we identify a narrow neighbourhood of the boundary dM with a part of dM x R_|_ 
and define a split of the cotangent bundle T*M = T*dM © T*(R). Let x = (x^), {i = 
1, 2, . . . , m — 1), be the local coordinates on dM and r be the normal geodesic distance 
to the boundary, so that N = dr = d/dr is the inward unit normal on dM. Near 
dM we choose the local coordinates x = (x^) = (x,r), (/i = 1,2, ...,m) and the split 
^ = (^^) = ((^uj) e T*M, where C = (0) e T*dM and to eR. With our notation, Greek 
indices run from 1 through m and lower case Latin indices run from 1 through m — 1. 

Our presentation differs from the one in since we always work with self-adjoint 
projectors. In this paper we are interested in the so-called generalized boundary con- 
ditions, when A is a first-order tangential differential operator acting on sections of the 
vector bundle Wq over dM. Any formally self-adjoint operator of first order satisfying 



the conditions (|2.24|) can be put in the form (hereafter, Vj denotes (m — l)-dimensional 



covariant differentiation tangentially, defined in ref. [^) 

A = (J - n) {^(rv, + v,r) + 5} (/ - n), (2.57) 

where F* G C°°{TdM ® E'n.d{Wo),dM) are some endomorphism- valued vector fields on 
dM, and S is some endomorphism of the vector bundle Wq, satisfying the conditions 

r = -r, S = S, (2.58) 



nF* = F^n = n5 = 5n = 0. 



(2.59) 
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Now we are going to determine under which conditions the boundary-value problem is 
strongly elliptic 

First of all, the leading symbol of the operator F should be elliptic in the interior of 
M. The leading symbol of the operator F reads 

= leP ■ / = 9^''{x)i,i, ■ J, (2.60) 

where ^ G T*M is a cotangent vector. Of course, for a positive-definite non-singular 
metric the leading symbol is non-degenerate for ^ 7^ 0. Moreover, for a complex A which 
does not lie on the positive real axis, A G C — R+ (R+ being the set of positive numbers), 

det (ai(F;x,0 - A) = (1^1' - A)'^'-^ ^ 0. (2.61) 

This equals zero only for ,^ = A = 0. Thus, the leading symbol of the operator F is elliptic. 

Second, the so-called strong ellipticity condition should be satisfied [|I], |I0|. As we 
already noted above, there is a natural grading in the vector bundles Wp and W'p which 
reflects simply the number of normal derivatives of a section of the bundle [|T|. The 
boundary operator Bp (p.25|) is said to have the graded order 0. Its graded leading symbol 
is defined by [g ^ 

.,(B,).(n (^_%,). (2.62) 

where, by virtue of (|2.57| ), 

T = -tapiA) = TXj, (2.63) 

( G T*dM being a cotangent vector on the boundary. By virtue of (|2.58|) the matrix T is 
anti-self-adjoint, 

f = -T, (2.64) 

and satisfies the conditions 

UT = TU = 0. (2.65) 

To define the strong ellipticity condition we take the leading symbol aL{F;x,r,(,^) of 
the operator F, substitute r = and u — > —idr and consider the following ordinary 
differential equation for a G C°°(V, dM x R+): 

K(F; x, 0, C, -tdr) - A ■ /] ipir) = 0, (2.66) 

with an asymptotic condition 

lim ip(r) = 0, (2.67) 

where ( G T*dM, A G C — R+ is a complex number which does not lie on the positive 
real axis, and {(, X) 7^ (0, 0). 

The boundary- value problem {F, B) is said to be strongly elliptic [|1^, p. 415] with 
respect to the cone C - R+ if for every C G T*dM, A G C - R+, ((, A) 7^ (0,0), and 
ip' G C°°(W',dM) there is a unique solution ip to the equation (|2.66| ) subject to the 
condition ( |2.67D and satisfying 



agiBp){xX)H^) = ^l^', 



(2.68) 
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with ilj{(p) G C°^{W,dM) being the boundary data defined by { ^T^ and (pi?D. 

A purely algebraic formulation of strong ellipticity can also be given, following Gilkey 
and Smith [1^. For this purpose, let us denote by Wp{C,\) the subsets of Wp corre- 
sponding to boundary data of solutions of Eq. ( |2.66| ) vanishing as r ^ ±00. Decompose 
(Tl(F; X, 0, C)i^) = Po<^'^ + Pi'^ + P2, where pj is homogeneous of order j in (. Then the 
differential equation ( f^.66| ) can be rewritten in the form of a first-order system 



— I 



dr + T{C,X) ( =0, (2.69) 



where 

r = zf JiM. (2.70) 

\Po (P2-A) Po PlJ 

Herefrom one sees that r does not have any purely-imaginary eigenvalues for {(, A) 7^ (0, 0) 
if {F, B) is strongly elliptic. It is then possible to re-express the strong ellipticity condition 
by saying that 

ag{BF){xX): W+{CA)^W'p (2.71) 

should be an isomorphism, for (C, A) ^ (0, 0), C e T*{dM), A G C - R+ [|IU|, p.416]. 
For a Laplace-type operator the equation (|2.66| ) takes the form 

" -5,2 + |Cr-A]y.(r) = 0, (2.72) 

where |CP = g^K^)CiCj- The general solution satisfying the decay condition at infinity, 
r ^ cx), reads 

V9(r) = xexp(-/ir), (2.73) 



where = y |CP — A. Since (C, A) 7^ (0,0) and A G C — R+, one can always choose 
Re/i > 0. 

The boundary data are now 

M^)=(4J- (2-74) 

Thus, the question of strong ellipticity for Laplace-type operators is reduced to the in- 
vertibility of the equations 



n \( X 



:2.75) 



for arbitrary Vo e C°°(iy^, 9M), ^'i e C'^{W[,dM). This is obviously equivalent to the 
algebraic criterion ( p.71|) . The eq. (|2.75|) can be transformed into (cf. 0) 

UX = ^'0, (2.76) 
(^/-iT)x = /i^o-^i. (2.77) 

Remember that ip'p = Pp'ipF with some ipp G C°°(VFi?, 9M). Therefore, (1 — n)?/^Q = 
Hip'i = and the first equation follows from the second one. Therefore, if the equation 
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( ^.77| ) has a unique solution for any ip^ and then the boundary- value problem is 
strongly elliptic. A necessary and sufficient condition to achieve this is expressed by the 
non-degeneracy of the matrix [fil — iT], i.e. 

det [ill - iT] ^ 0, (2.78) 



for any (C, A) ^ (0, 0) and A G C — R+. In this case the solution of eq. (|2.771 ) reads 



X={liI-iT)-^Wo-^[]. (2.79) 

Since the matrix iT is self-adjoint, it has only real eigenvalues, in other words the 
eigenvalues of are real and negative, < 0. It is clear that, for any non-real A G C — R, 



/i=Y|(^P — Ais complex and, therefore, the matrix [/iJ — iT] is non-degenerate. For real 
negative A, /i is real and we have /i > Thus, the condition ( p.78|) means that the 
matrix \C,\I — iT is positive-definite. 



A sufficient condition for that reads 



|C|/-iT>0. (2.80) 



|Cr/ + T'>0. (2.81) 



On the other hand, there holds, of course, \C,\'^I + T'^ < The eq. ( |2.81| ) means that 
the absolute values of all eigenvalues of the matrix {iT), both positive and negative, are 
smaller than whereas ( p.80| ) means that only the positive eigenvalues are smaller than 
Id, but says nothing about the negative ones. A similar inequality has been derived in 
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This proves the following theorem: 

Theorem 3 Let F be a Laplace-type operator defined by and Bp the generalized 

boundary operator given by ^2.2^ with the operator A defined by ^2.5% . Let C, G T*dM be 



a cotangent vector on the boundary and T = T^Q. The boundary-value problem (F,Bf) 
is strongly elliptic with respect to C — R+ if and only if for any non-vanishing C 7^ the 
matrix — iT is positive- definite, i.e. — iT > 0. 



2.3.2 Dirac-type operator 

The question of strong ellipticity for boundary-value problems involving Dirac-type oper- 
ators can also be studied. As is well known, the leading symbol of a Dirac-type operator 
reads 

ai(D;0 = -r(0. (2.82) 

Since the square of a Dirac-type operator is a Laplace-type operator, it is clear that the 
leading symbol aL{D;x,C,) is non-degenerate in the interior of M for ^ 7^ 0. Moreover, 
for a complex A one finds 

[aL{D; - A/]-' = MD; + A/] 777^- (2.83) 
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Therefore, [(Jl{D] ^) — XI] is non-degenerate when I^CP^-^^ 7^ 0. But this vanishes only for 
Im A = and Re A = ±|^| and arbitrary ^. Thus, for (^, A) ^ (0, 0) and A G C - R+ - R_, 
[aL{D;^) — A/] is non-degenerate. 

The boundary- value problem for a Dirac-type operator D is strongly elliptic if there 
exists a unique solution of the equation 

[zT{N)dr-T{C)-\I]^{r)=0, (2.84) 

for any (C, A) 7^ (0, 0), A ^ R+ U R__, such that 

lim ip(r) = 0, (2.85) 

r— >cxD 

PdMv) = ^'d, (2.86) 

for any -0^ e W'^^. 

The general solution of ( [2.84] ) satisfying the decay condition at infinity reads 

(/9(r) = xexp(-/ir), (2.87) 



where now fi = y |CP — A^. Again, since (C, A) 7^ (0, 0), A G C — R+ — R_, the root can 
be always defined by Re/i > 0. The constant prefactor x should satisfy the equation 

(X - \I)x = 0, (2.88) 

where 

X = aL{D- iiiN, C) = -iiiV{N) - r(C), (2.89) 
and the boundary condition 

PdX = ^'d, (2.90) 
with some ip'j^ = PdV^O; V^o ^ W/q- The eqs. ( p.88| ), ( |2.90| ) are reduced to 

X = ^'i, + X-, (2.91) 
[X + A(2P,, - /)] x^ = -{X- XI)^'j,. (2.92) 

Note that, since PdX- vanishes, the coefficient of Pd in Eq. ( p.92|) is arbitrary, and is 
set equal to 2 for convenience. Thus, the question of strong ellipticity for a Dirac-type 
operator is reduced to the matrix [X + A(2Pd — I)] being non-degenerate, 

det [X + X{2Pd - I)] ^ 0, (2.93) 

for any {(, A) 7^ (0, 0), A G C — R+ — R_. If this is satisfied, then the solution reads 

X = 2A [X + A(2P^ - /)]- V^. (2.94) 

Let us set (2Pd — I) = rj and let us compute the square of the matrix [X + Xrj]. For the 
boundary projectors defined by ( p.54| )-( p.56| ) the matrix t] is either rj = T{N), rj = r{u)/\u\ 
with some u G T*dM or (in even dimension m only) is exactly the chirality operator rj = C 
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(see ( |2.48D -( p.5(JD ). By using the properties ( ^.5[ ) of the Chfford algebra and those of the 
chirality operator rj (see (|2.51|) - (|2.53| )), we compute 



[X + \Cf = 2X^1, 



and 



X + \—T(u) 

\u\ 



- 2 

= 2A 






\u\ 



I. 



[2.95) 



(2.96) 



We see that, in both cases, the matrix [X + X{2Pd — I)] is non-degenerate for any {(, A) 7^ 
(0, 0) and any non- vanishing A that does not he on the real axis. 
Further, we find also 

[X + XT{N)f = 2\{\-tiJ,)I. (2.97) 



Bearing in mind that /i = y |CP "~ -^^j "we see that this does not vanish for any ( ^ 0. 
However, for C = and arbitrary A, with Im A > 0, this equals zero. Thus, the boundary 
projector Po = (/ ± r(A^))/2 does not lead to strong ellipticity, which only holds for 
ImA < 0. 

Thus, we have proven 

Theorem 4 The boundary-value problem {D, Bo), where D is a Dirac-type operator and 
Bd is a projector taking the forms i\2.5d{ ), or l \2. 5^ , or i\2. 5^ ), is strongly elliptic with 
respect to C — {0}, C — R and the lower half-plane Im A < 0, respectively. 



3 Resolvent and heat kernel 

Let w be a sufficiently large negative constant and A G C, Re A < w, be a complex number 
with a sufficiently large negative real part. Then the resolvent G{\) = {F — A/)^^ : 
L^(y,M) L'^{V,M) of the strongly elliptic boundary- value problem {F,B) is well 
defined. The resolvent kernel is a section of the tensor product of the vector bundles V 
and V* over the tensor-product manifold M x M, defined by the equation 

{F-\I)G{\\x,y) = 6{x,y) (3.1) 

with the boundary conditions 

BF^[Gi\\x,y)] = 0, (3.2) 

where S{x,y) is the covariant Dirac distribution, which is nothing but the kernel of the 
identity operator. Hereafter all differential operators as well as the boundary data map 
act on the first argument of the resolvent kernel (and the heat kernel), unless otherwise 
stated. This equation, together with the condition 

G{\\x,y) = G{\*\y,x), (3.3) 

which follows from the self-adjointness of the operator F, completely determine the resol- 
vent kernel. 
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Similarly, for t > the heat semi-group operator U{t) = exp{—tF) : L'^{V,M) 
L'^{V, M) is well defined. The kernel of this operator, called heat kernel, is defined by the 
equation 

{dt + F)Uit\x,y) = (3.4) 

with the initial condition 

U{0\x,y) = 6ix,y), (3.5) 

the boundary condition 

BFtlj[U{t\x,y)]=0. (3.6) 

and the self-adjointness condition 



U{t\x,y) = U{t\y,x). (3.7) 

As is well known [|I|, the heat kernel and the resolvent kernel are related by the Laplace 
transform: 

G{X\x,y)= J dte'^U{t\x,y), (3.8) 



w+ioo 

Uit\x,y) = ^^ J dXe~''GiX\x,y). (3.9) 

w—ioo 

Also, it is well known |]T| that the heat kernel U{t\x,y) is a smooth function near 

diagonal of M x M and has a well defined diagonal value U{t\x,x), and the functional 
trace 

Tt L2 exp{-tF) = I dYol{x)trvU{t\x,x). (3.10) 
Jm 

Moreover, the functional trace has an asymptotic expansion as t — 0+ 

Tr^2exp(-tF) ~ {Ant)-"^/' ^.t'^'^k^F, Bp), (3.11) 

fc>0 

and the corresponding expansion of the functional trace for the resolvent as A ^ — oo 
reads 



Ttl^ ( ^) G(A) ~ (47r)-™/2gr[(A;-m)/2+n+l](-A)(™-'=)/2-""iAfc/2(F,5F), (3.12) 



for n > 171/2. Here Ak/2{F, Bp) are the famous so-called {global) heat-kernel coefficients 
(sometimes called also Minakshisundaram-Plejel or Seeley coefficients). The zeroth-order 
coefficient is very well known: 

Ao = J dvo\ (x)tr yJ = vol (M) ■ dim{V). (3.13) 

M 

It is independent of the operator F and of the boundary conditions. The higher order 
coefficients have the following general form: 

Ak/2{F,Bp) = J dvol{x)tTvak/2{F\x)+ J dvo\{x)tT vbk/2{F, Bp\x), (3.14) 

M dM 
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where 0^/2 and hk/2 are the [local) interior and boundary heat-kernel coefficients. The 
interior coefficients do not depend on the boundary conditions Bp. It is well known 
that the interior coefficients of half-integer order, 0^+1/2, vanish The integer order 
coefficients are calculated for Laplace- type operators up to 04 [0. The boundary 
coefficients bk/2{F, Bp) do depend on both the operator F and the boundary operator Bp. 
They are far more complicated because in addition to the geometry of the manifold M 
they depend essentially on the geometry of the boundary DM. For Laplace-type operators 
they are known for the usual boundary conditions (Dirichlet, Neumann, or mixed version 
of them) up to 65/2 P, For generalized boundary conditions including tangential 
derivatives, to the best of our knowledge, they are not known at all. Only some special 



cases have been studied in the literature H, O, O] 



We are going to calculate below the next-to-leading coefficient Aii2{F, Bp) for the 
generalized boundary conditions. To do this, and also to study the role of the ellipticity 
condition, we will construct an approximation to the heat kernel U{t\x,y) near the di- 
agonal, i.e. for X close to y and for t — 0^. Since the heat kernel and resolvent kernel 
are connected by the Laplace transform, this is equivalent to studying an approximation 
to the resolvent kernel G{\\x,y) near the diagonal and for large negative A — —00 (this 
leads, in turn, to an approximate inverse of F — A/, called a parametrix) . 

Let us stress here that we are not going to provide a rigorous construction of the 
resolvent with all the estimates, which, for the boundary-value problem, is a task that 
would require a separate paper. For a complete and mathematically rigorous exposition 
the reader is referred to the classical papers ||10|, |T5|, [16], |17|, 0, [1^ . 

Here we keep instead to a pragmatic approach and will describe briefly how the ap- 
proximate resolvent kernel for A — > — 00, and hence the heat kernel for t — * 0+ can be 
constructed, and then will calculate both kernels in the leading approximation. This will 
allow us to compute the heat-kernel coefficient Ai/2. 

First of all, we decompose both kernels into two parts 

GiX\x,y)=G^i\\x,y) + GBiX\x,y), (3.15) 

U{t\x,y) = UUt\x,y) + UB{t\x,y). (3.16) 

Then we construct different approximations for Goo and Gb and, analogously, for Uoo and 
Ub- The ffist parts Goo(A|a;,j/) and Uoo{t\x,y) are approximated by the usual asymptotic 
expansion of the resolvent and the heat kernel in the case of compact manifolds without 
boundary when x ^ y, X —>■ —00 and t — > 0+. This means that effectively one introduces 
a small expansion parameter e reflecting the fact that the points x and y are close to 
each other, the parameter t is small and the parameter A is negative and large. This 
can be done by fixing a point x', choosing the normal coordinates at this point (with 
9fMu{x') = 5f,y) and scaling 



X — s> X 



' + e{x-x'), y ^ x' + e{y - x'), t ^ eH, X-^e~'^X (3.17) 



and expanding into an asymptotic series in e. If one uses the Fourier transform, then the 
corresponding momenta C, ^ T*M are large and scale according to 



(3.18) 
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This construction is standard []T| and we do not repeat it here. One can also use a 



completely covariant method [13, EO]. Probably the most convenient formula for the 



asymptotics as t 0"*", among many equivalent ones, is [|l^, 

U^{t\x,y) ~ {Ant)-^/'exp (-^) A^^ ^t^^a^, (3.19) 

^ fc>0 

where a = a{x,y) = r^{x,y)/2 is one half the geodesic distance between x and y, A = 
A{x,y) = g~^^'^{x)g~^^'^{y)det {—d^dya{x,y)) is the corresponding Van Vleck-Morette 
determinant, g = det gfj^^, and ak = ak{x,y) are the off-diagonal heat-kernel coefficients. 
These coefficients satisfy certain differential recursion relations which can be solved in 
form of a covariant Taylor series near diagonal |]13|. On the diagonal the asymptotic 
expansion of the heat kernel reads 

U^{t\x, x) ~ (47rt)-"^/2 ^ ^k^^^^^ ^3 20) 

fc>0 

The explicit formulas for the diagonal values of are known up to /c = 4 [jl3l. This 
asymptotic expansion can be integrated over the manifold M to get 

J dvo\{x)tivUooit\x,x) ~ {A-Kty^l^Y.^^ j dYo\{x)iivak{x,x). (3.21) 

Thus, integrating the diagonal of Uoo gives the interior terms in the heat-kernel asymp- 
totics. 

The asymptotic expansion of Goo(A|x, ?/) for A — > — oo is obtained from here by the 
Laplace transform ( |3.8| ) 

(2fc+2-m)/4 

M— J 

fc>0 



/ „ \ (2fc+2-m)/4 

Goo(A|a;, y) ~ {At,)-^'^ A"^ ^ 2 (^-^j irfe+i_„/2(v^=2A^)afc, (3.22) 



where K,,{z) = 7r/[2 sin(z/7r)][/_ 1,(2;) — 1,^(2;)] is the McDonald function (the modified Bessel 
function of third kind). It is singular on the diagonal x = y. However, for a sufficiently 
large n, n > 171/2, it becomes regular at the diagonal resulting in the asymptotic series as 
A —>■ —00 

y civol(x)try I J G'oo(A|x,?/) 

~ (47r)-™/2 ^r(A;-- + n+l) {-X)f~k-n~^ f ^^^i ^^-^^r v ak{x, x). (3.23) 



x=y 



k>0 



For a strongly elliptic boundary-value problem the diagonal of the boundary part 
UB{t\x,x) is exponentially small, i.e. of order ~ exp(— r^(x)/t), where r{x) is the normal 
geodesic distance to the boundary, as t 0"*" if x ^ dM. So, it does not contribute 
to the asymptotic expansion of the heat-kernel diagonal outside the boundary as t ^ 
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0+. Therefore, the asymptotic expansion of the total heat-kernel diagonal outside the 
boundary is determined only by Uoo 

U{t\x, x) ~ (47rt)-™/2 J2 i^(^k{x, x), x^ dM. (3.24) 

A:>0 

The point is, the coefficients of the asymptotic expansion of the diagonal of the boundary 
part UB{t\x,x) as t — > 0^ behave near the boundary like the one-dimensional Dirac 
distribution 5{r{x)) and its derivatives. Thus, the integral over the manifold M of the 
boundary part UB{t\x,x) has an asymptotic expansion as t — 0+ with non-vanishing 
coefficients in form of the integrals over the boundary. Therefore, it determines the local 
boundary contributions hk/2 to the heat-kernel coefficients Ak/2- It is well known that the 
coefficient y4i/2 is a purely boundary contribution [Q. It is almost obvious that it can 
be evaluated by integrating the fibre trace of the boundary contribution Ub oi the heat 
kernel in the leading order. 

Of course, this approximation is obtained without taking into account the boundary 
conditions. Therefore, Goo satisfies approximately the equation (|3.1|) but does not satisfy 



the boundary conditions ( |3.2| ). This means that the compensating term GB{^\x,y) is 
defined by the equation 

(F-A/)Gb(A|x,i/) = (3.25) 

with the boundary condition 

Bpi^ [G'oo(A|x, y) + Gb{\\x, y)] = 0. (3.26) 
Analogously, UB{t\x,y) is defined by 

{dt + F)UB{t\x,y) = (3.27) 

with the initial condition 

UB{0\x,y) = 0, (3.28) 

and the boundary condition 

BFtlj[U^it\x,y) + UBit\x,y)] = 0. (3.29) 

The most difficult problem is to find the compensating terms Gb(A|x, y) and [/s(t|x, y). 
These functions are important only near the boundary where they behave like distributions 
when t — > 0"*" or A — *■ —oo. Since the points x and y are close to the boundary the 
coordinates r{x) and r{y) are small separately, hence not only the difference [r{x) —r{y)] 
but also the sum [r{x) + r{y)] is small. This means that we must additionally scale 
r(x) — * er{x), r{y) Br{y). By contrast, the point x' is kept fixed on the boundary, so 
the coordinates x' do not scale at all: x' x'. 

Thus, we shall scale the coordinates x = {x,r{x)) and y = {y,r{y)), the parameters t 
and A and momenta ^ = (C, uj) G T*M, with ( G T*dM and u G T*R, according to 



X 



x+e{x~x), y—^x+e{y — x), r(x) — er(x), r{y) er{y), (3.30) 
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—A, C-^-C, u; ^ - 

E'^ e e 

The corresponding differential operators are scaled by 



t^eh, A -.-A, C-^tC, ^^-y- (3-31) 



d -(9, dr -<9r, dt \du d\ e^dx- (3.32) 

We will call this transformation just scaling and denote the scaled objects by an index e, 
e.g. G^. The scaling parameter e will be considered as a small parameter in the theory 
and we will use it to expand everything in power series in e. We will not take care about 
the convergence properties of these expansions and take them as formal power series. In 
fact, they are asymptotic expansions as e — >■ 0. At the very end of calculations we set 
e = 1. The non-scaled objects, i.e. those with e = 1, will not have the index e, e.g. 
G'|j|e=i = Gb- Another way of doing this is by saying that we will expand all quantities 
in the homogeneous functions of the coordinates {x — ct'), {y — x'),r{x),r{y), the momenta 
^ = {(, Lj) and the parameters t, X. 

First of all, we expand the scaled operator in power series in e 

F'r^Y. ^^~'Fn, (3.33) 

n>0 

where F„ are second-order differential operators with homogeneous symbols. The bound- 
ary operator requires a more careful handling. Since half of the boundary data (|2.8|) 
contain normal derivatives, formally ipo = ^\r=o and ipi = drip\r=o, (^.9|), would be of 
different order in e. To make them of the same order we have to assume an additional 
factor e in all ipi G C°°{Wi,dM). Thus, we define the graded scaling of the boundary 
data map by 

^•(-)-(.frJ,)-(4l:*)i:!J^^'-'- 

so that the boundary data map ip does not scale at all. This leads to an additional factor 
e in the operator A determining the boundary operator Bp ( p.25| ). Thus, we define the 
graded scaling of the boundary operator by 

^^=(.A* /-V)- t^-^^) 
which has the following asymptotic expansion in e: 

B^pr^Y.^'^Bp^n), (3.36) 

n>0 

where Bp(^n) are first-order tangential operators with homogeneous symbols. Of course, 

Fo = - d\ (3.37) 

BFm = (^l /-Ho)' ^^-^^^ 
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where 

= g^\x')djdk, Ao = T\x')dj, Uo = U{x'). (3.39) 

Note that all leading-order operators Fq, Bf{o) and Aq have constant coefficients and, 
therefore, are very easy to handle. This procedure is called sometimes "freezing the 
coefficients of the differential operator" . 

The subsequent strategy is rather simple. Expand the scaled resolvent kernel in e 



G'L~E^'""''"Goo(„), (3.40) 

n>0 
n>0 

and substitute into the scaled version of the equation ( |3.25| ) and the boundary condition 
( p.26| ). Then, by equating the like powers in e one gets an infinite set of recursive equations 
which determine all GB{n)- The zeroth-order term Gb(o) is defined by 

(Fo - A/)Gb(o) = 0, (3.42) 

and the boundary conditions, 

Bpi^)^ (^00(0) + G'b(o)) = 0, (3.43) 
where Fq and -B_p(o) are defined by ( |3.37| )-( |3.39| ). The higher orders are determined from 

k 

(Fo - XI)GBik) = - E FnGBik-n). = 1, 2, . . . , (3.44) 

n=l 

and the corresponding boundary conditions, 

k 

Bf{0)'^ (Goo{k) + GB{k)) = — E BF{n)i' (G'oo(fe-n) + G B{k-n)) ■ (3.45) 

n=l 

The Goo(n) are obtained simply by expanding the scaled version of ( |3.22p in power series 
in e. The operator Fq is a partial differential operator, but fortunately, has constant 
coefficients. By using the Fourier transform in the boundary coordinates {x — x') it 
reduces to an ordinary differential operator of second order. This enables one to find 
easily its resolvent. We will do this in next subsection below. Before doing this, let us 
stress that the same procedure can be applied to get the boundary part UB{t\x,y) of the 
heat kernel. 



3.1 Leading order resolvent 

In this subsection we determine the resolvent to leading order, i.e. Go = Goo(o) + Gb(q)- 
As we already outlined above, we fix a point x' G DM on the boundary and the normal 
coordinates at this point (with gik{x') = 6ik), take the tangent space TdM and identify 
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the manifold M with Mq = TdM x R+. By using the exphcit form of the zeroth-order 
operators Fq, Bq and Aq given by ( p.37D -( p73^ ) we obtain the equation 

(-dl -d^-\) G,{\\x,y) = 5{x - y), (3.46) 

and the boundary conditions 

noGo(A|x,i/)L^^^^^ = 0, (3.47) 

(/ - Ho) {dr + iVidj) Go(A|x, y) [^^^^^ = 0, (3.48) 

where Ho = n(x'), Fq = T^{x'). Hereafter the differential operators always act on the first 
argument of a kernel. Moreover, for simplicity of notation, we will denote IIo and Fq just 
by n and F-' and omit the dependence of all geometric objects on x' . To leading order 
this cannot cause any misunderstanding. Further to this, the resolvent kernel should be 
bounded, 

lim G,{\\x,y)= lim G'o(A|x, y) = 0, (3.49) 

and self- adjoint, 

GoiX\x,y) = Go{X*\y,x). (3.50) 

Since the above boundary-value problem has constant coefficients, by using the Fourier 
transform in {x — y) 

G(o)(A|x,y)= j ^|-^e<-(^-*)G(o)(A|C,r(x),r(t/)), (3.51) 

where C, € T*dM, ( ■ x = CjX\ we obtain an ordinary differential equation 

{-d', + \Cf- a) Go(A|C,r(x),r(y)) = 6[p{x,y)], (3.52) 
where |CP = g'^{x')CiCj, 

p{x,y) = r{x) - r{y), (3.53) 

the boundary conditions 

nGo(A|C,0,r(|/)) = 0, (3.54) 
(/ - n) (dr + iT) Go(A|C, r{x)Ay)) , , n = 0' (3.55) 

r(x)=\j 

and self-adjointness condition 



Go(A|C,r(x),r(y)) = Go(A*|C, r(y), r(x)). (3.56) 
Here T is an anti-self-adjoint matrix 

T = T(C) = F ■ C = F^O, f = -T, (3.57) 

satisfying the conditions 

(/ - H)T = T{I - H) = T. (3.58) 
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The first 'free' part of the solution of this problem Goo(o), which gives the resolvent 
kernel on the line R, is almost obvious. Since the equation (p.52|) has constant coefficients 
and there are no boundary conditions, its solution depends only on the difference p(x, y) = 
r{x) — r{y) 



2p 
1 

2^ 



J_g-M|p(a;,2/)|_ 



where 

/i = \/lCp-A, Re/x>0, 
and 6 is the usual step function 



(3.59) 
(3.60) 



9{x) = 



1 if X > 
if X < 0. 



(3.61) 



The boundary part Gb{o) takes into account the boundary conditions ( p.54| ) and ( p.55| ). 
On requiring regularity at infinity and self-adjointness of the kernel, we can write down 
the solution up to a constant self-adjoint matrix h = h, 



2/i 



(3.62) 



/^From eqs. ( |3.59| )- (|3.62|) we find 

G,{X\C,0,r{y)) = ^{I + h)e-^^-^y\ 

ZjJ, 



drGo{X\C,r{x),r{y)) 



r(x)=0 



-{I-h)e 



-tJ.riy) 

(3.63) 

Taking into account the conditions ( |3.57| )-( P755| ) we reduce the boundary conditions to 
the form 

n(/ + /i) = 0, (3.64) 

{I -U)fi{I -h)+iT{I + h) = 0. (3.65) 

Since h is self-adjoint and Til = IIT = it must have the form h = all + (/ — n)/5(T)(J — 
n), where a is a constant and /? is a matrix depending on T. Substituting this into ( |3.64| ) 
we immediately obtain a = — 1. From the second boundary condition (|3.65| ) we get 



(fxl -iT)p = pi + iT. 



(3.66) 



Since the problem is assumed to be strongly elliptic, the matrix pi — iT is invertible for 
any {(, A) 7^ (0, 0). Thus, we obtain the solution of the boundary conditions in the form 



hip;C) = -n + {I -n){pl -tT)-\pI + tT){I -U) 
= I -2U + 2iT{pI -iT)-\ 



(3.67) 
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In a particular case of mixed Dirichlet and Neumann boundary conditions, when P = 0, 
this reduces to 

h = I~2n. (3.68) 

Note that the function h{fi, () depends actually on the ratio (/ fi,m other words, it satisfies 
the homogeneity relation 

= /i(/x,C). (3.69) 

More generally, 

Go(^^A|-^C,Vtr(x),Vtr(|/)^ = v^Go(A|C,r(x),r(y)) (3.70) 

This holds for Goo(o) as well as for Gb{o)- As a function of A, G(o) is a meromorphic 
function with a cut along the positive real axis from to oo. The "free" part 6*00(0) has 
no other singularities, whereas the boundary part Gb{o) has simple poles at the eigenvalues 
of the matrix + T^. Remember that the strong ellipticity condition |2.80| requires these 



eigenvalues to be positive and smaller than i.e. < |CP + < \(\'^. 
3.2 Leading order heat kernel 

Using Go{X\x,y), we can also get the zeroth-order heat kernel Uo{t\x,y) by the inverse 
Laplace transform 

w+ioo 

Uo{t\x,y) = ^ / dXe~'^GoiX\x,y) (3.71) 



dC ""T dX 



(27r)'"-i J 2TTi 

1 ii) — ioo 



e 



-*W(-y)Go(A|C,r(x),r(y)) (3.72) 



where w is a negative constant. Now, by scaling the integration variables X ^ X/t and 
C C/Vt and shifting the contour of integration over X {w ^ w/t, which can be done 
because the integrand is analytic in the left half-plane of A) and using the homogeneity 
property we obtain immediately 



ui+ioo 



R™-1 w—ioo 

(3.73) 

Next, let us change the variable A according to 

A = |Cr + ^', (3.74) 

where u = ifi, and, hence, Imci; > 0. In the upper half-plane, Imuj > 0, this change 
of variables is single valued and well defined. Under this change the cut in the complex 
plane A along the positive real axis from to 00, i.e. ImA = 0, |CP < ReA < 00 
is mapped onto the whole real axis Ima; = 0, —00 < Reu; < -|-oo. The interval 
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ImA = 0, < ReA < ICp on the real axis of A is mapped onto an interval Keu = 
0, < Imuj < Id, on the positive imaginary axis of u. As a function of u the resolvent 
Go is a meromorphic function in the upper half plane, Imct; > 0, with simple poles 
on the interval Retu = 0, < Imu < \(\, on the imaginary axis. The contour of 
integration in the complex plane of is a hyperbola going from (e*^'^/^)oo through the 



point Lu = -y/lCP -w to (e*^/^)oo. It can be deformed to a contour C that comes from 
—oo + ie, encircles the point u = i\(\ in the clockwise direction and goes to +oo + ie, 
where e is an infinitesimal positive parameter. The contour C does not cross the interval 
Rea; = 0, < Imu; < |C|, on the imaginary axis and is above all the singularities of the 
resolvent. 

After such a tranformation we obtain 

Uoit\x,y) = (Ant)-/^ [ :J^e-\<\'^^^<''yy^' 



m — l 



X /-.-,(-,.)«„(|cr,.=|C,llf,f ) (3.75) 
Substituting here Go we find 

X 

c 

The first "free" part is easily obtained by computing the Gaussian integrals over u and 
(. We get 

Uo{t\x, y) = (47rt)— /2 exp f-i |x - y\'') + f/B(o) y), (3.77) 



4t' 

where 

f/^(o)(t|x,y) = (4^t)-/^ / ^^(^e-l^l^+<-(^-*)/^ (3.78) 



/■ doj 

X 



/ ^ exp |-.^ + ,,M:^^)±!M| {(J _ 2n) - 2T{. I + rr} . 

Here, the part proportional to (/ — 211) also contains only Gaussian integrals and is easily 
calculated. Thus (cf. pT[| ) 

UBio){t\x, y) = (Ant)-^/' jexp - y\' - ^[r(a:) + r{y)Yj (/ - 2n) + n{t\x, y)} , 

(3.79) 

where 

xr-C(u;/ + r-C)^\ (3.80) 
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where we substituted T = T ■ ( = T^Q. 

Herefrom we obtain easily the diagonal value of the heat kernel 



U^o){t\x,x) = {4nt)-^/'{l + 



exp 



(/-2n) + <l>(r(a;)/Vt)|, 



f3.811 



where 



-2 / r.c(./ + r.c)-'. (3.82) 

R™-i C ^ 



Now, by using the representation (remember that Im (cj / + F ■ C) > 0) 
and computing a Gaussian integral over u, we obtain 



^{z) = 2 
Integrating by parts over p we get 



^(m-l)/2 



^pe-ICP-(p+^)^_^p2ipr.c 



dp 



(3.83) 



(3.84) 



oz 



where 



^iz) 



dC 



R" 



7r(™-i)/2 Jq 



dpe 



-ICP-(p+2)2+2ipr.c 



(3.85) 



(3.86) 



It is not difficult to show that, as 2; ^ 00, the functions '^{z) and $(z) are exponentially 
small: 



^ ' 2z 



$(2) ~ 

z^ 



-r^ + o(2 



-2^ 



where = (yfjjTT-'. For 2; = 0, by using the change C — "we obtain 



*(0)^f / 



^(m-l)/2 



-ICP-(r-C)2 



(3.87) 
(3.88) 

(3.89) 



R" 



Note that this integral converges when the strong ellipticity condition (F ■ (^)^ + |Cp > 
is satisfied. 

Now, we take the diagonal f/(o)(t|a;, x) given by ( 3.81 ), and integrate over the manifold 
M. Because the boundary part \Jb{q) is exponentially small as rix) ^ cxd we can in fact 
integrate it only over a narrow strip near the boundary, when < r(x) < 5. The difference 
is asymptotically small as t — 0"*". Doing the change of variables z = r/\/t we reduce the 
integration to < z < 6/ \/t. We see, that as t ^ 0"*" we can integrate over z from to 00. 
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The error is asymptotically small as t — > 0+ and does not contribute to the asymptotic 
expansion of the trace of the heat kernel. 
Thus, we obtain 

Tr L2Uoit) ~ (47rt)-'"/2 + ^^^^ ^ . . (g gg) 

where Aq is given by (|3.13|) and 

Ai/2= J dvo\{x)tTvai/2, (3.91) 

dM 

with 



«i/2 = ^(/-2n + A/2), (3.92) 
2 °° 

f3y2 = ^ f dz<l>{z). (3.93) 





Now, using (|3.85|) and ( |3.89|) and the fact that ^(oo) = we get easily 



Pi/2 = -2 + ^^(0) 

dc 



^(m-l)/2 



W ^e-^^%C^ - l) , (3.94) 



where 
Eventually, 



0-1/2 



/ - 2n + 2 



Further calculations of general nature, without knowing the algebraic properties of the 
matrices A^^, seem to be impossible. One can, however, evaluate the integral in form of 
an expansion in the matrices A^^, or F*. Using the Gaussian integrals 

,-\c\\ c --^or a- . (3 97) 

we obtain 

A/2 = 2 Y: i-^r7^^9^.^2 " ' ' 9^2.-.^2J^^' " " " r^^"\ (3-98) 

n>i y^-) ^ 

and, therefore. 



«i/2 = ^{^ - 2n + 2|: (-1)" (^l^^n. ■ ■ ■ ^.„-..„r(- . . . F-) 



(3.99) 



Since our main result ( p.94|) is rather complicated, we now consider a number of 
particular cases of physical relevance. 
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1. The simplest case is, of course, when the matrices A^^ vanish, A^^ 
gets the famihar result for mixed boundary conditions Q 



Pi/2 = 0, ai/2 



(/-2n). 



0. One then 



(3.100) 



2. The first non-trivial case is when the matrices F* form an Abelian algebra, 

[T\P] = 0. 

One can then easily compute the integral ( |3.94| ) explicitly. 



/5l/2 



Therefore 



«l/2 



(/ + r2)-V2_i 



{/ - 2n + 2 [(/ + r2)-i/2 - /] I 



(3.101) 



(3.102) 



(3.103) 



In the case 11 = 0, this coincides with the result of ref. 0, where the authors consid- 
ered the particular case of commuting F* matrices (without noting this explicitly). 

3. Let us assume that the matrices A^'' = F'^-^T'^^ form an Abelian algebra, i.e. 

[A^^A'"]=0. (3.104) 



The evaluation of the resulting integral over ( yields 



1/2 



\detTaM[S', + A'j])-'/'-I 



(3.105) 



where the determinant det tsm is taken over the indices in the tangent space to the 
boundary (we used det^ = !)• By virtue of ( p. 92] ) one gets the final result 



0-1/2 



I -2U + 2 



{detTdM[S'j+A\ 



-1/2 



(3.106) 



4. A particular realization of the last case is when the matrices F* satisfy the Dirac-type 
condition 



so that 



FT-'' + Pr = 2g 



im — 11 



:9 



*ip2 



(3.107) 



(3.108) 



(m — 1) 

Then, of course, the matrices A^^ commute and the result is given by eq. (|3.106|) . 
But the determinant is now 



m—l 



det TdM[S'j + A\ 



(m — 1 ) 



(3.109) 
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Thus, 



ai/2 



/ - 2n + 2 



-(■m-l)/2 



/ + 



(3.110) 



Note that this differs essentially from the resuh of ref. and shows again that 
the resuh of ref. |^ apphes actually only to the completely Abelian case, when all 
matrices P commute. Note also that, in the most interesting applications (e.g. in 
quantum gravity), the matrices F* do not commute |]T2[. The result ( |3.94| ), however, 
is valid in the most general case. 

5. A very important case is when the operator A is a natural operator on the boundary. 
Since it is of first order it can be only the generalized Dirac operator. In this case 
the matrices P satisfy the condition 



rP + Pr = 2A'i = -2 Kg'^{i - n). 



(3.111) 



where k is a constant. Hence the matrices A^^ obviously commute and we have the 
case considered above (see ( p.l06| )). The determinant is easily calculated. 



detTaAf[5} + 4] =n+(J-n)(l 
and we eventually obtain 



\in—l 



0-1/2 



^{-n + (/-n) 



2(1 



^(m-l)/2 



(3.112) 



(3.113) 



Thus, a singularity is found at k = 1. This happens because, for k = 1, the strong 
ellipticity condition is violated (see also [|lll). Indeed, the strong ellipticity condition 



T' + icr/ = (r ■ 0' + icr/ = icnn + (1 -«:)(/ - n)] > 0, 

implies in this case k < 1 (cf. [|ll 



(3.114) 



4 Analysis of ellipticity in a general gauge theory on 
manifolds with boundary 

In this section we are going to study gauge-invariant boundary conditions in a general 
gauge theory (for a review, see |^). A gauge theory is defined by two vector bundles, 

V and G, such that diml^ > dimG. V is the bundle of gauge fields ip G C°°{V, M), and 
G is the bundle of parameters of gauge transformations e G C°°{G, M). Both bundles 

V and G are equipped with some Hermitian positive-definite metrics E, E'^ = E, and 7, 
7"'' = 7, and with the corresponding natural scalar products (, )y and (, )g- 

The infinitesimal gauge transformations 



6ip = Re 



(4.1) 
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are determined by a first-order diff'erential operator R, 

R : C°^{G,M) ^C°^{V,M). (4.2) 
Further, one introduces two auxiliary operators, 

X : C^iV, M) C°°(G', M) (4.3) 

and 

F : C~(G',M) ^ C~(G',M), (4.4) 
and one defines two difi^erential operators, 

L = XR:C'^{G,M)^C'^{G,M) (4.5) 

and 

H = XYX: C^{V,M)-^C^{V,M), (4.6) 

where X = E~^X'^'y. The operators X and Y should satisfy the following conditions (but 
are otherwise arbitrary): 

1) The differential operators L and H have the same order. 

2) The operators L and H are formally self-adjoint (or anti-self-adjoint). 

3) The operators L and Y are eUiptic. 

From these conditions we find that there are two essentially different cases: 

Case I. X is of first order and Y is of zeroth order, i.e. 

X = R, Y = Ig, (4.7) 
where R = ^~^R^E. Then, of course, L and H are both seconc?- orc?er differential operators, 

L^RR, RR. (4.8) 

Case II. X is of zeroth order and Y is of first order. Let TZ be the bundle of maps of 
G into V, and let /3 e 7?. be a zeroth-order differential operator. Then 

X^$, Y^^R, (4.9) 

where (3 = ^~'^f3^E, and the operators L and H are of first order, 

L^PR, (3PRP = (3Lp. (4.10) 

We assume that, by suitable choice of the parameters, the second-order operator RR can 
be made of Laplace type and the first-order operator fiR can be made of Dirac type, and, 
therefore, have non-degenerate leading symbols, 

det GTLiRR) 0, (4.11) 
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detGai(/?i?) ^0. (4.12) 

The dynamics of gauge fields ip G C°°(K M) at the hnearized (one- loop) level is described 
by a formally self-adjoint (or anti-self-adjoint) differential operator, 

A : C^{y, M) C~(V, M). (4.13) 

This operator is of second order for bosonic fields and of first order for fermionic fields. 
In both cases it satisfies the identities 

Ai? = 0, i?A = 0, (4.14) 

and, therefore, is degenerate. 

We consider only the case when the gauge generators are linearly independent. This 
means that the equation 

OLiRY = 0, (4.15) 
for ^ 7^ 0, has the only solution e = 0. In other words, 

Ker(7i(i?) = 0, (4.16) 

i.e. the rank of the leading symbol of the operator it! equals the dimension of the bundle 



rank(TL(-R) = dimG. 


(4.17) 


We also assume that the leading symbols of the generators R 


are complete in that they 


generate all zero-modes of the leading symbol of the operator 


A, i.e. all solutions of the 


equation 




(7i(A)^ = 0, 


(4.18) 


for ^ 7^ 0, have the form 




if = (TL{R)e, 


(4.19) 


for some e. In other words. 




Kerai(A) = {aL(i?)e | e G G}, 


(4.20) 


and hence 




rank crj^(A) = dim V — dim G. 


(4.21) 


Further, let us take the operator H of the same order as the operator A and construct a 


formally (anti-)self-adjoint operator. 




F = A + i/. 


(4.22) 


so that 




a^(F)=az.(A) + a^(i/). 


(4.23) 


It is easy to derive the following result: 





Proposition 1 The leading symbol of the operator F is non- degenerate, i.e. 

detyai(F;x,0 7^0, (4.24) 

for any ^ 7^ 0. 
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Proof. Indeed, suppose there exists a zero-mode, say (po, of the leading symbol of the 
operator F, i.e. 

aL{F)ipo = ^oaUF) = 0, (4.25) 

where (f = Lp'^E. Then we have 

^0<yL{F)aL{R) = v?o^L(Xr)ai(L) = 0, (4.26) 
and, since (Jl{L) is non-degenerate, 

ifioaLiXY) = aL{YX)^o = 0. (4.27) 

But this implies 

aLiH)ifo = 0, (4.28) 

and hence 

aLiF)ifo = aL{A)ipo = 0. (4.29) 

Thus, ipo is a zero-mode of the leading symbol of the operator A, and according to the 
completeness of the generators R must have the form ipo = ai^R)^ for some e. Substituting 



this form into the eq. ( [4.27|) we obtain 

aL{YX)aL{R)e = aL{YL)e = 0. (4.30) 

Herefrom, by taking into account the non-singularity of aiiYL), it follows e = (po = 0, 
and hence the leading symbol of the operator F does not have any zero-modes, i.e. it is 
non-degenerate. 

Thus, the operators L and F have, both, non-degenerate leading symbols. In quantum 
field theory the operator X is called the gauge-fixing operator, F the gauge-field operator, 
the operator L the (Faddeev- Popov) ghost operator and the operator Y in the Case II the 
third (or Nielsen-Kallosh) ghost operator. The most convenient and the most important 
case is when, by suitable choice of the parameters it turns out to be possible to make both 
the operators F and L either of Laplace type or of Dirac type. The one-loop effective 
action for gauge fields is given by the functional superdeterminants of the gauge-field 



operator F and the ghost operators L and Y |22 



r = hog (Sdet F) - log (Sdet L) - hog (Sdet Y). (4.31) 
4.1 Bosonic gauge fields 

Let us consider first the case of bosonic fields, when A is a second-order formally self- 
adjoint operator. The gauge invariance identity ( [4.14|) means, in particular. 



aL{A)aLiR) = 0. (4.32) 

Now we assume that both the operators L = RR and F = A + RR are of Laplace type, 
i.e. 

aURR) = lei'/c, (4.33) 
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aL{F) = aLiA) + aLiRR) = \^\'L 



(4.34) 



On manifolds with boundary one has to impose some boundary conditions to make these 
operators self-adjoint and elliptic. They read 



BL^ie) = 0, 



(4.35) 



BfHv) = 0, (4.36) 

where '?/'(e) and are the boundary data for the bundles G and V, respectively, and 
B]^ and Bp are the corresponding boundary operators (see section 2.1). In gauge theories 
one tries to choose the boundary operators Bl and Bp in a. gauge- invariant way, so that 
the condition 

Bp^jjiRe) = (4.37) 

is satisfied identically for any e subject to the boundary conditions ([4.35|) . This means 
that the boundary operators B^ and Bp satisfy the identity 



Bp[i;,R]{I-BL)=0, 



(4.38) 



where [ip, R] is the commutator of the linear boundary data map ip and the operator R. 

We will see that this requirement fixes completely the form of the as yet unknown 
boundary operator Bp- Indeed, the most natural way to satisfy the condition of gauge 
invariance is as follows. Let us decompose the cotangent bundle T*M in such a way that 
C, = {N, () G T*M, where is the inward pointing unit normal to the boundary and 
( e T*dM is a cotangent vector on the boundary. Consider the restriction Wq of the 
vector bundle V to the boundary. Let us define restrictions of the leading symbols of the 
operators R and A to the boundary, i.e. 



U = aLiA;N) 

u = aL{R; N) 

/i = apiR; C) 
From eq. ([4.32|) we have thus the identity 

Uu = 0, 

Moreover, from (f4.33 ) and ( 4.34| ) we have also 

uu = Ig, 
ufi + fiu = 0, 

fif^ = ICI%, 
U = I -uu. 



9M' 
dM' 



(4.39) 

(4.40) 
(4.41) 

(4.42) 

(4.43) 

(4.44) 
(4.45) 
(4.46) 
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From ( |4.42| ) and ( |4.43| ) we find that 11 : Wq — >• Wq is a self-adjoint projector orthogonal 
to u, 

= n, Uu = o, n = n. (4.47) 



Then, a part of the boundary conditions for the operator F reads 



dM 



The gauge transformation of this equation is 



URe 



dM 



0. 



0. 



(4.48) 



(4.49) 



The normal derivative does not contribute to this equation, and, therefore, if Dirichlet 
boundary conditions are imposed on e. 



dM 



0, 



(4.50) 



the equation ( |4.49| ) is satisfied identically. 

The easiest way to get the other part of the boundary conditions is just to set 



dM 



0. 



(4.51) 



Bearing in mind eq. ( |4.5| ) we find that, under the gauge transformations 
transformed into 



Le 



dM 



0. 



Tl), this is 



(4.52) 



If some e is a zero-mode of the operator L, i.e. e G Ker (L), this is identically zero. For all 
e ^ Ker (L) this is identically zero for the Dirichlet boundary conditions ([4.50|) . In other 
words, the requirement of gauge invariance of the boundary conditions ( ^4.36| ) determines 
in an almost unique way (up to zero- modes) that the ghost boundary operator Bl should 
be of Dirichlet type. Anyway, the Dirichlet boundary conditions for the operator L are 
sufficient to achieve gauge invariance of the boundary conditions for the operator F. 

Since the operator R in the boundary conditions ( |4.51| ) is a first-order operator, the set 
of boundary conditions (|4.48| ) and ( [4.51|) is equivalent to the general scheme formulated 
in section 2.1. Separating the normal derivative in the operator R we find exactly the 
generalized boundary conditions ( p.l5| ) with the boundary operator Bp oi the form ( |2.25| ) 
with a first-order operator A : C°^{WQ,dM) C^^iW^^dM), the matrices F-' being of 
the form 

F^' = -w^h. (4.53) 



These matrices are anti-self-adjoint, F* = — F*, and satisfy the relations 

Yir = ru = 0. 



(4.54) 



Thus, one can now define the matrix 



(4.55) 
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where /x = /i-^'Cj, and study the condition of strong elhpticity ( p.80|) . The condition of 
strong elhpticity now reads 

\C\I-iT=\C\I + iJ^ufiu>0. (4.56) 
Further, using the eqs. ([4.53|) , ([4.44]) and ( [4.45| ) we evaluate 



A'^ = r^T^) = -(/ - n)/iVV - n). (4.57) 

Therefore, 

= A'XiCj = -{I - n)/x/x(/ - n), (4.58) 

and 

+ \cfi = ICI'n + (/ - n)[|c|'/ - M(/ - n). (4.59) 

Since for non-vanishing ( the part proportional to 11 is positive-definite, the condition of 
strong elhpticity for bosonic gauge theory means 

(/-n)[|cr/-/x/x](/-n) >o. (4.60) 

We have thus proved a theorem: 

Theorem 5 Let V and G be two vector bundles over a compact Riemannian manifold 
M with smooth boundary, such that dimV > dimG. Consider a bosonic gauge theory 
and let the first-order differential operator R : C°°{G, M) — > C'^{V,M) be the gener- 
ator of infinitesimal gauge transformations. Let A : C°°(V,M) C°^(V,M) be the 
gauge-invariant second-order differential operator of the linearized field equations. Let the 
operators L = RR : C'^{G,M) G'^{G,M) and F = A + RR be of Laplace type 
and normalized by (Tl{.L) = I^^/g- Let ai^R-^N) = v and (7l{R-X) = fi be the 



dM 



dM 



restriction of the leading symbol of the operator R to the boundary, N being the normal 
to the boundary and ( G T*dM being a cotangent vector, and U = I — uu. 

Then the generalized boundary-value problem {F, Bp) with the boundary operator Bp 
determined by the boundary conditions ( \4.43i ) and l \4 . 51\ ) is gauge-invariant provided that 
the ghost boundary operator B^ takes the Dirichlet form. Moreover, it is strongly elliptic 
with respect to C — R+ if and only if the matrix [\C\I + ivv^v] is positive- definite. A 
sufficient condition for that reads 

(j-n)[|cp/-M/i](/-n)>o. (4.61) 



4.2 Fermionic gauge fields 

In the case of fermionic gauge fields, A is a first-order formally self-adjoint (or anti-self- 
adjoint), degenerate (i.e. gauge- invariant) operator with leading symbol satisfying 

aL{A)aL{R) = 0. (4.62) 

Now we have the case II, hence, we assume both the operators L = /3R and F = A+(3/3R/3 
to be of Dirac type, i.e. the operators and are of Laplace type: 



(4.63) 
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az.(A)+/5^az.(/?)^] =1^1'/. (4.64) 

Note that now we have two systems of Dirac matrices, 7^ on the bundle G and on V. 
They are defined by the leading symbols of the Dirac-type operators L and F, 

a^iL) = -ri,. ^l{F) = -V^^,. (4.65) 

Let us define 



b = aLiA;N) 



<9M' 



c = aL{A;C) 



divr 



(4.66) 



where C, € T*dM. Bearing in mind the notation ([4.40|) and ( [4.41|) , we have from ( [4.62| ) 



bu = 0, cu + bfi = 0. 
Moreover, from (|4.63|) we have also 

Pufju = Ig, 
j3vl3n + /5/i/5z/ = 0. 



Similarly, from (|4.64|) we get 

{b + (3^v^f 

(b + /5^z//3) (c + + (c + (3^^^) (b + (3^u^) = 0. 

Herefrom, we find that the following operators: 



Pv = \[I-T{N)] = ^{I + b + l3^u^), 



(4.67) 

(4.68) 
(4.69) 

(4.70) 
(4.71) 

(4.72) 
(4.73) 

(4.74) 



are projectors, 

= Pg, P^ = Pv 

The boundary conditions for the Dirac-type operators L and F are given by projectors 
(see (pi])). 



Pr.e 



dM 



Ppip 



dM 



0, 



0. 



(4.75) 
(4.76) 



The problem is to make them gauge- invariant. Here the projectors Pl and Pp are defined 
by means of the matrices 7^ and F^, respectively. The gauge transformation of the 
equation (|4.76| ) is 

Noting that, on the boundary. 



dM 



0. 



(4.77) 



Re 



dM 



<9A/' 



(4.78) 



and assuming that yU commutes with P^, we get herefrom two conditions on the boundary 
projectors, 

Pfi^MIg - Pl) = 0, (4.79) 

Pf{I - uPuP){Ig - Pl) = 0. (4.80) 

Such gauge-invariant boundary operators always exist. We will construct them explicitly 
in section 6 in the course of studying the Rarita-Schwinger system. 
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5 Strong ellipticity in Yang— Mills theory 

The first physical application that wc study is the strong eUipticity condition in Yang- 
Mills theory. Now G = A is the Lie algebra of a semi-simple and compact gauge group, 
and V is the bundle of 1-forms taking values in A, i.e. V = T*M (g) G. Let h be the 
Cartan metric on the Lie algebra defined by 

hab = —C'^ad.C'^cb: (5.1) 

C"bc being the structure constants of the gauge group, and the fibre metric E on the 
bundle V be defined by 

E{ip, if) = -tTAgi^p, ip), (5.2) 

or, in components, 

E^a'b = g^'hab. (5.3) 

The Cartan metric is non-degenerate and positive-definite. Therefore, the fibre metric 
E is always non-singular and positive-definite. Henceforth we will suppress the group 
indices. The generator of gauge transformations is now just the covariant derivative 

R^V^, R^-tTgV^. (5.4) 

The leading symbols of these operators are 

(tl{R; = i^iG, <jl{R; = -i^i, (5.5) 

where ^ = tr is a map ^ : T*M R. First of all, we see that 

<7l(L; i) = aL{RR-^ i) - I^Ig = \i?lG. (5.6) 

ULiH- = ULiRR-, = e ® (5.7) 

so that the operator L — RR is indeed a Laplace-type operator. The gauge-invariant 
operator A in linearized Yang-Mills theory is defined by the leading symbol 

aL{/^) = {\i\^-i®i)l. (5.8) 

Thus, the operator F = A -|- if is of Laplace type, 

aL{F)=aL{^ + H) = \i\^I. (5.9) 

Further, we find 

1/ = i NIc, u = -itr gN, (5.10) 

and 

Ii = iC,Ig, fj, = -ttrgC (5.11) 

where N is the normal cotangent vector, and C, G T*dM. 
The projector 11 has the form 

n = (5.12) 
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where 

q = l-N(g)N. (5.13) 

In components, this reads 

q^-^^5;-N,N\ (5.14) 
Thus, the gauge-invariant boundary conditions are 



dM 



0, (5.15) 



g'-'V,ip.^^^ = 0. (5.16) 
Since (N = NC =< (, N >= 0, we find from (plO|) and ( pTID 

ftu = ufi = 0, (5-17) 

and hence the matrices F* in (|2.57|) , as well as T = TXi, vanish: 

p = 0, T = 0. (5.18) 

Therefore, the matrix + C^-^ = ICP-^ is positive-definite, so that the strong ellipticity 
condition (|2.80|) or ([4.61| ) is satisfied. 
Thus, we have 

Theorem 6 A Laplace-type operator F acting on 1-forms taking values in a semi-simple 
Lie algebra G, F : C°°(T*M®G, M) C'^{T*M®G, M), with the boundary conditions 
(\5.13[ ) and ( \5.1di) , is elliptic. 

Since such boundary conditions automatically appear in the gauge-invariant formula- 
tion of the boundary conditions in one-loop Yang-Mills theory, we have herefrom 

Corollary 1 The boundary-value problem in one-loop Euclidean Yang-Mills theory de- 
termined by a Laplace-type operator and the gauge-invariant boundary conditions defined 
by ( |5. i5[ ) and l \5.1(^ , is strongly elliptic with respect to C — R_|_. 



6 Ellipticity for the Rarita— Schwinger system 

The next step is the analysis of Rarita-Schwinger fields. The bundle G is now the bundle 
of spinor fields taking values in some semi-simple Lie algebra A, i.e. G = S ® A, and V is 
the bundle of spin-vector fields, in other words, V is the bundle of 1-forms taking values 
in the fibre of G, i.e. V = T*M ® G. The whole theory does not depend on the presence 
of the algebra A, so we will omit completely the group indices. 

Let 7^ be the Dirac matrices and 7 be the Hermitian metric on the spinor bundle S 
determined by 



YY + Yr = 2g'"'Is, 7^ = -£7^ 7^ = 7, 



(6.1) 
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where 7^ = 7 ^7^"''^. The fibre metric E on the bundle V is defined by 



where 



E{ip,ip) = ^l^E'^''ip,, (6.2) 



with a parameter a. By using ( |6.1|) it is easily seen that E is Hermitian, i.e. 

E^"" = E''^, (6.4) 

if a is real. The inverse metric reads 

7-1-1 

e;. = 9,. - (TT^^'^^- (6.5) 

Therefore, the fibre metric E is positive-definite only for a > —1/m and is singular for 
a = —1/m. Thus, hereafter a ^ —1/m. In fact, for a = —1/m the matrix Efj_'^ becomes 
a projector on a subspace of spin-vectors (p^ satisfying the condition 7^v^^ = 0. 

The generator of gauge transformations is now again, as in the Yang-Mills case, just 
the covariant derivative 

R = 6V^, (6.6) 
where 6 is a normalization constant, with leading symbol 

aL{R;0=ib^lG. (6.7) 

Now we have the Case II of section 4, and hence we define the map jS : G ^ V and its 
adjoint /3 : — G by 

= tE-,^^''^ = 1 = 77V., (6.8) 
^ 6(1 + am) 

so that 

_ — STTl 

0^(1 -|- am) 



The operator X of Eq. (|4.3| ) is now equal to P so that the operator L = PR is of Dirac 
type with leading symbol 

aL{L) = -r^^. (6.10) 
The leading symbol of the operator H = P/3R/3 reads 

aL{H)^x = jTTT^ ,lxrY^,Vu. (6.11) 

o"^(l + am) 

The gauge-invariant operator A is now the Rarita-Schwinger operator with leading symbol 

aL{A)^x = eE^^Y^''^,^,. (6.12) 
Here and below we denote the antisymmetrized products of 7-matrices by 

^Ml-Mn ^ . . .^Mn]_ (6.13) 
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Of course, the leading symbol is self-adjoint and gauge- invariant, in that 



<7l(A) = ai,(A), (7i,(A)(7i(i?) = 0, (6.14) 

where crL(A) = E~^al^{A)E. Further, the leading symbol of the operator F — A + H 
reads 

Er^-f'^f^" + 7x7^7'' 



Using the properties of the Clifford algebra we compute 



CuVu- (6.15) 



[ [l + ma) V b^J [l + ma) J 

(6.16) 

Moreover, one can prove the following property of representation theory: 



Proposition 2 The matrices r^{0) defined by 

{T^mx'' = r^x%0) ^ + -7A7V + (^ - - (^ + 1)7A^^^ (6.17) 

TTli 

where cu is defined by 

= , 6.18 

m 



form a representation of the Clifford algebra, i.e. 



ry (o)rv(o) + ^V(o)^^''(o) = 2gi'^5'is. (6.i9) 



It is thus clear that the set of matrices 

t\f{a') = r-^''(a')fV(0)T/(a'), (6-20) 

with arbitrary non-degenerate matrix T{a') depending on a parameter a', also forms a 
representation of the Chfford algebra. By choosing 

T/(a')=5f + «W, (6.21) 

and, hence, 

T-\^{a') = 5i - ,, 7A7^ (6.22) 
(1 -I- ma') 

we prove, more generally, what follows. 
Corollary 2 The matrices T^{a') defined by 

{f^{a'))x'' = f\^{a') ^ 5^7^ + ci7A7^7'^ + C26^xY + Cs^xg'", (6.23) 

where 

ci = a' H , 6.24 

(l + ma') V 2 ; ' ^ ^ 
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C2 = -(m -2-mu) (a' + ^^^j , (6.25) 

2 / , uj + l\ , , 

(1 + ma') V 2 / 

wzt/i a' an arbitrary constant a' ^ —1/m, and uj defined by ( 6.18J , form a representation 
of the Clifford algebra, i.e. satisfy the eq. ^6.1^\ ). 

Note that, by choosing a' = — (c<j + l)/2, all the constants Ci, C2 and C3 vanish. 

The operator F with leading symbol ( |6.16| ) should be of Dirac type. Thus, by imposing 
that the term in curly brackets in Eq. (|6.16|) should coincide with the right-hand side of 
Eq. (|6.23|) , one finds a system whose solution is 



a 



[m - 4 + (m - 2)uj] 



m — 4 



a 



^Jm-2 

The simplest case is m = 4, one has then a = a' = 0, 6 = ±^/2. 

After this choice the operator F is of Dirac type, i.e. cr^lF) = —eV^^^, with 



(6.27) 

(6.28) 
(6.29) 



[m 



(m 



:lxg 



(6.30) 



Thus, we have two Dirac-type operators, L and F, which have elliptic leading symbols. By 
choosing the appropriate boundary conditions with the projectors Pl and Pp the system 
becomes elliptic. The problem is to define the boundary projectors in a gauge-invariant 



way. 



Let Pl the boundary projector for the ghost operator L, 



Pi.e 



dM 



0. 



(6.31) 



Remember that it satisfies the symmetry condition (|2.44| ). Then we choose the boundary 
conditions for the gauge field in the form 



dM 



dM 



0, 



where g is a projector defined by ( |5.13| ). The gauge transformation of ( |6.32| ) reads 



PLq/^ne 



dM 



(6.32) 
(6.33) 

(6.34) 



and does not include the normal derivative. We are assuming that the projector P^ com- 
mutes with the tangential derivative (as is usually the case, we find that their commutator 
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vanishes identically by virtue of the boundary conditions on e). The gauge transformation 
of eq. ( |6.33D is proportional exactly to the operator L, 



dM 



PLLe 



dM 



0. 



(6.35) 



By expanding e in the eigenmodes of the operator L we find that this is proportional to the 
boundary conditions ( p. 31 ) on e, and therefore vanishes. Thus, the boundary conditions 
( |6.32| ) and ( |6.33|) are gauge-invariant. They can be re- written in another convenient form, 



This defines eventually the boundary projector Pp for the gauge operator F, 
Pf/ = [5/ - N.N^^N'' - N^N^YY]PL[qu^ + N,l% 



(6.36) 



(6.37) 



If the projector Pl satisfies the symmetry condition ( p.44| ) then so does the projector Pp 
(of course, one has to check it with the matrix F'^A*'^). 
Thus, we have shown that 

Theorem 7 The boundary-value problem for the RaritaSchwinger system with the bound- 
ary conditions ( |6'. 31\ )- ^6. is gauge-invariant and strongly elliptic provided that the pro- 
jector Pl satisfies the condition ( ^.44 )- Particular examples of such projectors are given 

by gM)-(\iM)- 



7 Euclidean quantum gravity 

Generalized boundary conditions similar to the ones studied so far occur naturally in 
Euclidean quantum gravity [§, ^. The vector bundle G is now the bundle of cotangent 
vectors, G = T*M, and V is the vector bundle of symmetric rank-2 tensors (also called 
symmetric 2-forms) over M: V = T*M\/T*M, V being the symmetrized tensor product. 
The metric on the bundle G is, naturally, just the metric on M, and the fibre metric E 
on the bundle V is defined by the equation 

where a is a real parameter. One has also, for the inverse metric, 

E ^ab cd = ga{cgd)b — TTT- zQabgcd- (7.2) 

(1 + am) 

We do not fix the a parameter from the beginning, but rather are going to study the 
dependence of the heat kernel on it. It is not difficult to show that the eigenvalues of the 
matrix E are 1 (with degeneracy m(m + l)/2 — 1) and (1 + am). Therefore, this metric 
is positive-definite only for 

a > - — , (7.3) 

m 
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and becomes singular for a = —1/m. Thus, hereafter a ^ —1/m. 

The generator R of infinitesimal gauge transformations is now defined to be the Lie 
derivative of the tensor field along the vector field e, 



The adjoint generator R with the metric E is defined by its action, 
The leading symbol of the ghost operator L = RR takes now the form 



aLiL; Oea = 2E,"^%^,e, = 5'J^\' + (1 + 2^)^' 



eft. 



(7.4) 



(7.5) 



(7.6) 



Therefore, we see that it becomes of Laplace type only for a = —1/2. Note also that the 
operator L has positive-definite leading symbol only for a > — 1, and becomes degenerate 
for a = —1. Further, the leading symbol of the operator H = RR reads 



-.cd 



^cd- 



(7.7) 



The gauge- invariant operator A is well known (see, e.g. |2^, ^). It has the following 
leading symbol: 

Thus, we see that the operator F = A + is of Laplace type only in the case a = —1/2. 
Let us, however, consider for the time being a Laplace-type operator E on symmetric 
rank-2 tensors with a fibre metric ( [7.1| ) with an arbitrary parameter a. 
Further, we define the projector 11 



n = Uab"'' = Qia" qbf - 



a 



a + 1 



cd 



(7.9) 



where qab = Qab — NaNb. It is not difficult to check that it is self-adjoint with respect to 
the metric E, i.e. 11 = 11, and that 



tryll = —m{m — 1). 
2 



(7.10) 



Thus, we consider a Laplace-type operator E acting on symmetric rank-2 tensors with 
the following boundary conditions: 



^ab'^'^'^cd 



dM 



0, 



cd 



dM 



0. 



(7.11) 
(7.12) 
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Separating the normal derivative we find from here the boundary operator Bp oi the 
form (|2.25|) , the operator A being given by ( |2.57|) with the matrices defined by 



= r 



cd 



1 



ab 



-—NaN,e'^'N^^ + Nf^^elN'N''. (7.13) 
1 + a ' 



It is not difficult to check that these matrices are anti-self-adjoint and satisfy the conditions 
and (^391) . The matrix T = T ■ C reads (cf. g) 

T = -7T^^Pi+P2, (7.14) 
(1 + a) 

where 

= pi,,^-^ = N^N,C^'^N^\ P2 = P2ab"' = N^aCb)N''N'', (7.15) 

and (a = e^Ci; so that CaN"' = 0. It is important to note that 

nT = rn = o. (7.16) 

We also define further projectors, 

p = p^,-^ ^ ^^N^^C,^N^X''\ (7.17) 

p = p,,''' = NaN,N'^N'', (7.18) 

which are mutually orthogonal: 

pp = pp = 0. (7.19) 

The matrices pi and p2, however, are nilpotent: ^1=^2 = 0' ^^"^ their products are 
proportional to the projectors: 

P1P2 = l\C\'p, P2P1 = ^ICIV (7.20) 



Therefore, we have 



where 



We compute further 



Last, by using 



we obtain 



= -7^r^,\(\"(P + P) = -^"(P + P)' (7-21) 
2(1 + a) 

I / J CI- (7.22) 

^2(1 + a) 



T2" = (2r)2"(p + p), (7.23) 

T2"+i = (ir)2"T. (7.24) 

trp = trp = l, (7.25) 

trpi = trp2 = 0, (7.26) 



trT^" =2(zr)2", tr T2"+^ = tr T = 0. (7.27) 
This suffices to prove: 
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Lemma 1 For any function f analytic in the region \z\ < r, one has 



f{T) 



/(O) [/ - p - p] + ^[/(zr) + fi-tT)]{p + p) 
+^[f{^r)-f{-^r)]T, 



tr/(T) 



m{m + 1) 



/(O) + f{^T) + f{-iT) 



(7.28) 
(7.29) 



Thus, the eigenvalues of the matrix T are 



spec (T) 



with degeneracy 
ir with degeneracy 1 
—IT with degeneracy 1 



m(m+l) 
2 



(7.30) 



with T defined by eq. ( \7.23i ). 



This means that the eigenvalues of the matrix for a non-vanishing Q are and 
— 1/[2(1 + a)]|Cp. Thus, the strong elhpticity condition ( p. 81 ), which means that the 
matrix (T^ + (^I)^ for non-zero should be positive-definite, takes the form 



1 



2(1 + a) 



+ 1 > 0. 



(7.31) 



This proves eventually 



Theorem 8 The boundary-value problem for a Laplace-type operator acting on sections 
of the bundle of symmetric rank-2 tensors with the boundary conditions l \7.11\ ) and l \7.1^ ) 
is strongly elliptic with respect to C — R_|_ only for 



a > —- 



(7.32) 



Remarks. First, let us note that the condition ( [7.32 ) of strong ellipticity is compatible 
with the condition ( |7.3| ) of positivity of the fibre metric E. Second, it is exactly the 
value a = —1/2 that appears in the gauge- invariant boundary conditions in one-loop 
quantum gravity in the minimal DeWitt gauge. For a general value a ^ —1/2, the 
operator F resulting from A and H is not of Laplace type, which complicates the analysis 
significantly. In other words, we have a corollary H: 



Corollary 3 The boundary-value problem in one-loop Euclidean quantum gravity, with a 
Laplace-type operator F acting on rank-two symmetric tensor fields, and with the gauge- 



invariant boundary conditions ( \7.1]\ ) and (\7.12J , the fibre metric being E with a = —1/2, 
is not strongly elliptic with respect to C — R+. 
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We can also evaluate the coefficient ai/2 of heat-kernel asymptotics. This is most easily 
obtained by using the representation (|3.94|) of the coefficient /3i/2 in form of a Gaussian 



integral. Using eq. ( |7.21| ) we get ffist 



exp(-r^ 



exp 



2(1 + a 



(7.33) 



Therefore, from eq. (|3.94|) we obtain 



1/2 



R" 



7r(™-i)/2 



exp 



1 + 2a , 2 



2(1 + a) 



(p + p). 



(7.34) 



One should bear in mind that p is a projector that depends on ( (see ( [7.17| )). Although 
p is singular at the point C = 0, the integral is well defined because of the difference of 
two exponential functions. Calculating the Gaussian integrals we obtain 



13: 



1/2 



'2{a + iy 
1 + 2a 



(m-l)/2 



P + 



{m — 1) 



TC 



ai/2 = ^ <; / - 2n + 2 



'2{a + l) 
l + 2a 



(m-l)/2 



P + 



m — 1 



where is yet another projector: 



(7.35) 



(7.36) 



(7.37) 



Last, from eq. ( p.92|) , by using the traces of the projectors (|7.1CI| ), (|7.25| ) and (|7.37|) , and 
the dimension of the bundle of symmetric rank-2 tensors, 



we obtain 



dim V = tvyl = -m{m + 1), 
2 



/TT I 1 

tr v'ai/2 = < --m{m - 3) - 4 + 4 



"2(a + !)■ 




_ 1 + 2a _ 


(m-i)/2| 



(7.38) 



(7.39) 



Thus, we see that there is a singularity at a = —1/2, which reflects the lack of strong 
ellipticity in this case. 



7.1 Heat-kernel diagonal in the non-elliptic case 

Consider now the case a = —1/2. From eq. ( [7. 281 ) we have, in particular. 



(/p-^T)-i 



-(/-p-p) 



p p^ 
-(/ - p - p) - y(p + P) - i^T. 

jjL A 



^ -Xp + p)+^ ^ 



X' 



(7.40) 
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For the boundary-value problem to be strongly elliptic, this determinant should be non- 
vanishing for any ((^, A) 7^ (0, 0) and A G C — R+, including the case A = 0, C 7^ 0. 
Actually we see that, for any non-zero A ^ R+, this determinant does not vanish for 
any C,. However, for A = and any C, it equals zero, which means that the corresponding 
boundary conditions do not fix a unique solution of the eigenvalue equation for the leading 
symbol, subject to a decay condition at infinity. This is reflected by the simple fact that 
the coefficient ai/2 of the asymptotic expansion of the heat kernel is not well defined, 
in that the integrals that determine it are divergent. At the technical level, the non- 
ellipticity is reflected in the fact that the heat-kernel diagonal, although well defined, has 
a non-standard non-integrable behaviour near the boundary, i.e. for r — > 0. 

To prove this property, let us calculate the fibre trace of the heat-kernel diagonal. 
From ( 3.81| ) we have 



where 



tTvUo{t\x,x) = (47rt)-™/2 + cie"'''/* + J{r/Vi) 

mim + 1) 



Co 



dim 



ci = iiv{I - 2n) 



m(m — 3) 



J{z) 



dC 



7r(™-i)/2 



duj 



-\(\'^-uj^+2iujz 



c 



tivT-aui+T-cr 



(7.42) 

(7.43) 
(7.44) 
(7.45) 



Now, for a = —1/2 the parameter r ( |7.22|) determining the eigenvalues of the matrix 
T = T ■ ( is equal to r = |^|. By using ( [7.29| ) we get 



Jiz) = -4 



dC 



R" 



7r(™-i)/2 



du 



C 



TC 



(7.46) 



Remember that the contour C lies in the upper half plane: it comes from —00 + ie, 
encircles the point u = in the clockwise direction and goes to 00 + ie. The integral 
over uj is calculated by using the formula 



00 

J dtu f{uj) = —27TiRes^=i\(\ f{uj) + J dcu /{to). 



(7.47) 



The integrals over ( can be reduced to Gaussian integrals by lifting the denominator in 
the exponent (cf. (|3.83|) ) or by using spherical coordinates. In this way we get eventually 



J{z) = 2(m - l)r(m/2) - 2(m - 1) / rfuu^/^-^e""' 

2/ J 



(7.48) 



This can also be written in the form 



J{z) = 2{m-l)z'"' r(m/2) -7(m/2,z2) = 2(m - l)z-"r(m/2, 2^) 



(7.49) 
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by using the incomplete 7-functions 



00 



7(a,x) =y ciuu"-V", T{a,x)=J duu''~^e~''. (7.50) 
^ 

It is immediately seen that the function J{z) is singular as 2; — > 

J(^)~2(m-l)r(m/2)-^. (7.51) 

By using the identity 

r(a, X) = x'^-^e"^ + (a - l)T{a - 1, x) (7.52) 
we also find that, when z — ^ 00, J{z) is exponentially small, i.e. 

J{z) r^2{m-l)z-^e-'\ (7.53) 

Note that, for a = — |, one has tryT^ = — 2|C|^, and hence tryFT-' = —2g'^^ and 
tr^r^ = — 2(m — 1). Thus, we see that the asymptotics ( |7.53| ) as 2; ^ 00 corresponds to 
eq. dpI). 

The singularity at the point z = results exactly from the pole at u = In the 
strongly elliptic case all poles lie on the positive imaginary line with Imu < so that 
there is a finite gap between the pole located at the point with the largest value of the 
imaginary part and the point 

Thus, we obtain 



tivUo{t\x,x) = {Ant)-""/^ 



Co + cie-'-'/* + 2(m - 1) l-^j T{m/2, r^/t) 



(7.54) 



Here the first term is the standard first term in the heat-kernel asymptotics which gives the 
familiar interior contribution when integrated over a compact manifold. The boundary 
terms in ( [7.54| ) are exponentially small as t — 0^, if r is kept fixed. Actually, these 
terms should behave, as t — > 0"*", as distributions near the boundary, so that they give 
well defined non-vanishing contributions (in form of integrals over the boundary) when 
integrated with a smooth function. The third term, however, gives rise to an unusual 
singularity at the boundary as r ^ (on fixing t): 

tr vUo{t\x, x) (47r)-'"/22(m - l)r(m/2)-^. (7.55) 

This limit is non-standard in that: i) it does not depend on t and ii) it is not integrable 
over r near the boundary, as r ^ 0. This is a direct consequence of the lack of strong 
ellipticity. 
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8 Concluding remarks 

We have studied a generalized boundary-value problem for operators of Laplace type, 
where a part of the field is subject to Dirichlet boundary conditions, and the remaining 
part is subject to conditions which generalize the Robin case by the inclusion of tangential 
derivatives. The corresponding boundary operator can be always expressed in the form 
( [^.25| ), where A is a tangential differential operator of the form ( |2.57[ ). The fermionic 
analysis for Dirac-type operators has also been developed. The strong ellipticity of the 
resulting boundary-value problem is crucial, in particular, to ensure the existence of the 
asymptotic expansions frequently studied in the theory of heat kernels [|1|. In physical 
problems, this means that the one-loop semi-classical expansions of the Green functions 
and of the effective action in quantum field theory are well defined and can be computed 
explicitly on compact Riemannian manifolds with smooth boundary. The occurrence of 
boundaries plays indeed a key role in the path-integral approach to quantum gravity , 



and it appears desirable to study the strong ellipticity problem for all gauge theories of 
physical interest, now that a unified scheme for the derivation of BRST-invariant boundary 
conditions is available 0. 

We have thus tried to understand whether the following requirements are compatible: 

(i) The gauge-field operator, F, should be of Laplace type (and the same for the ghost 
operator); 

(ii) Local nature of the boundary operator Bp (in fact, we have studied the case when Bp 
is a first-order differential operator, and boundary projectors for fermionic fields); 

(iii) Gauge invariance of the boundary conditions; 

(iv) Strong ellipticity of the boundary- value problem (F, Bp). 

First, we have found a condition of strong ellipticity (see (|2.80| )) for a generalized boundary- 
value problem for a Laplace-type operator. For operators of Dirac type, one finds instead 
the strong ellipticity properties described in our Theorem 4. Second, we have constructed 
the resolvent kernel and the heat kernel in the leading approximation (sect. 3.1 and sect. 
3.2) and evaluated the first non-trivial heat-kernel coefficient A1/2 in Eqs. (|3.96|) and 
( ^.99| ). Third, we have found a criterion of strong ellipticity of a general gauge theory in 
terms of the gauge generators (see ( [4.61| )). As physical applications of the above results 
we have studied the Yang-Mills, Rarita-Schwinger and Einstein field theories. Interest- 
ingly, only in the latter the strong ellipticity condition is not satisfied if the conditions 
(i), (ii) and (iii) hold. Moreover, the gauge- invariant boundary conditions for the Rarita- 
Schwinger system have been found to involve only the boundary projector for the ghost 
operator. As far as we know, our results as well as the consequent analysis of the phys- 
ical gauge models, are completely original, or extend significantly previous work in the 
literature |, |, O]. 

Since we find that, for gravitation, the four conditions listed above are not, in general, 
compatible, it seems that one should investigate in detail at least one of the following 
alternatives: 
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(1) Quantum gravity and quantum supergravity on manifolds with boundary with 
gauge-field operators which are non-minimal; 

(2) Non-local boundary conditions for gravitation and spin-3/2 fields 



(3) Non-gauge-invariant 'regularization' of the boundary conditions, which suppresses 
the tangential derivatives and improves the ellipticity; 

(4) Boundary conditions which are not completely gauge-invariant, and hence avoid the 
occurrence of tangential derivatives in the boundary operator. 

The latter possibility has been investigated in [^] and has been widely used in the 
physical literature (see [^] and references therein). The first 3 options are not yet (com- 
pletely) exploited in the literature. In particular, it is unclear whether one has to resort 
to non-minimal operators to preserve strong ellipticity of the boundary- value problem. 
Moreover, non-local boundary conditions which are completely gauge- invariant, compat- 
ible with local supersymmetry transformations at the boundary (cf. [^), and ensure 
strong ellipticity, are not easily obtained (if at all admissible). Indeed, it should not be 
especially surprising that gauge theories have, in general, an essentially non-local char- 
acter. Since the projectors on the physical gauge-invariant subspace of the configuration 
space are non-local, a consistent formulation of gauge theories on Riemannian manifolds 
(even without boundary) is necessarily non-local (see p3|). 

Thus, there is increasing evidence in favor of boundaries raising deep and unavoidable 



foundational issues for the understanding of modern quantum field theories I^Uj. The 
solution of such problems might in turn shed new light on spectral geometry and on the 
general theory of elliptic operators [IT], pO|. 
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